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Z)*±(2317) and KD scattering from B° s decay 
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We study the B° —> DJ (DK) + weak decay, and look at the DK invariant mass distribution with 
the aim of obtaining relevant information on the nature of the D*q (2317) resonance. We make a 
simulation of the experiment using the actual mass of the D *q" (2317) resonance and recent lattice 
QCD relevant parameters of the KD scattering amplitude. We then solve the inverse problem 
of obtaining the KD amplitude from these synthetic data, to which we have added a 5% or 10% 
error. We prove that one can obtain from these ’’data” the existence of a bound KD state, the 
KD scattering length and effective range, and most importantly, the KD probability in the wave 
function of the bound state obtained, which was found to be largely dominant from the lattice QCD 
results. This means that one can obtain information on the nature of the D* q" (2317) resonance from 
the implementation of this experiment, in the line of finding the structure of resonances, which is 
one of the main aims in hadron spectroscopy. 

PACS numbers: 


I. INTRODUCTION 

The very narrow charmed scalar meson D* 1 g (2317) was 
first observed in the D+ tt° channel by the BABAR Col¬ 
laboration 1] and its existence was confirmed by CLEO 
Q, BELLE [3[ and FOCUS [|[ Collaborations. Its mass 
was commonly measured as 2317 MeV, which is approxi¬ 
mately 160 MeV below the prediction of the very success¬ 
ful quark model for the charmed mesons Q. Due to its 
low mass, the structure of the meson D*±( 2317) has been 
extensively debated. It has been interpreted as a cs state 
6 HJ- two-meson molecular state jTl|-[l8j . K—D- mixing 
13 |, four-quark states iiHIl or a mixture between two- 
meson and four-quark states (24[. Additional support to 
the molecular interpretation came recently from lattice 
QCD simulations [25l - l28i ]. In previous lattice studies of 
the D* 0 (2317), it was treated as a conventional quark- 
antiquark state and no states with the correct mass (be¬ 
low the KD threshold) were found. In Refs. [25|,[27j, with 
the introduction of KD meson operators and using the ef¬ 
fective range formula, a bound state is obtained about 40 
MeV below the KD threshold. The fact that the bound 
state appears with the KD interpolator may be inter¬ 
preted as a possible KD molecular structure, but more 
precise statements cannot be done. In Ref. [26j] lattice 
QCD results for the KD scattering length are extrapo¬ 
lated to physical pion masses by means of unitarized chi¬ 
ral perturbation theory, and by means of the Weinberg 
compositeness condition [29|, |3Cj the amount of KD con¬ 
tent in the D* 0 (2317) is determined, resulting in a sizable 
fraction of the order of 70% within errors. A reanalyis of 
the lattice spectra of Refs. [25j,[27} has been recently done 


in Ref. [Hf, going beyond the effective range approxima¬ 
tion and making use of the three levels of Refs. 

As a consequence, one can be more quantitative about 
the nature of the D s o(2317), which appears with a KD 
component of about 70%, within errors. 

In addition to these lattice results, and more precise 
ones that should be available in the future, it is very 
important to have some experimental data that could be 
used to test the internal structure of this exotic state. 

Here we propose to use the experimental KD invariant 
mass distribution of the weak decay of ► Dj {DK) +1 
in order to obtain information about the internal struc¬ 
ture of the _D*(f (2317) state. There are not yet exper¬ 
imental data for the decay — > D~{DK) + . How¬ 

ever, since the branching fractions for the decays B° —> 
D* + D*~ and B% -A Df D*~ + D* + Dj are respectively 
1.85% and 1.28%, we believe that the branching frac¬ 
tion for the —t DjD^ decay, should not be so dif¬ 
ferent from that and it will be seen through the chan¬ 
nel — > D~ (DK) + . This is why it is really impor¬ 

tant to have theoretical predictions for the DK invari¬ 
ant mass distribution that considers the formation of 
the D*q"( 2317) state. At this point, it is worth stress¬ 
ing that recently, in the reactions B° —> D~D°K + and 
B + -A D°D°K + studied by the BABAR Collaboration 
[3lj . an enhancement in the invariant DK mass in the 
range 2.35 — 2.50 GeV is observed, which could be asso¬ 
ciated with this 11*1^(2317) state. It is also interesting 
to quote that in a different reaction, —> D 0 K~n + 1 

the LHCb Collaboration also finds an enhancement close 
to the KD threshold in the D°K~ invariant mass dis¬ 
tribution, which is partly associated to the D* 0 (2317) 
resonance (Hj. 
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1 Throughout this work, the notation (DK) + refers to the isoscalar 

combination D°K + + D+K°. 
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FIG. 1: Mechanism for the decay B° —> D g (DK) + . 


In Fig. □ we show the mechanism for the decay B° s —> 
Dj{DK) + . The idea is to take the dominant mechanism 
for the weak decay of the B® into Dj plus a primary cs 
pair. The hadronization of the initial cs pair is achieved 
by inserting a qq pair with the quantum numbers of the 
vacuum: uu + dd + ss + cc, as shown in Fig. |T] There¬ 
fore, the cs pair is hadronized into a pair of pseudoscalar 
mesons. This pair of pseudoscalar mesons is then allowed 
to interact to produce the D*g (2317) resonance, which is 
considered here as mainly a DK molecule [HJ. The idea 
is similar to the one used in Ref. |33j for the formation of 
the / o (980) and / o (500) scalar resonances in the decays 
of B° and B° s . 

The paper is organized as follows. In Section [II] we 
settle the formalism for our study. Namely, in Subsec¬ 
tion QTA] we study the ( DK) + elastic scattering ampli¬ 
tude, and in Subsection IIIBI we study the differential 
decay width for the process B ° —y D~ {DK) + . As said 
before, there is not yet experimental information con¬ 
cerning the differential decay width for this process. For 
this reason, we will have to generate synthetic data for 
this decay in order to explore if this reaction is suitable 
for the study of the ( DK) + final state interactions and 
the D*(/(2317) bound state. The generation and analy¬ 
sis of these synthetic data, which constitutes the results 
of the work, are done in Subsection m Conclusions are 
delivered in Section HVl 


II. FORMALISM 


In this work the influence of the presence of the 
D*o'(2317) in the process B / —> D~{DK) + is investi¬ 
gated. The H*g'(2317) is considered mainly as a bound 
state of the DK system, so we address the elastic DK 
scattering amplitude in Subsection III Al Then, the differ¬ 
ential decay width for the B° s —> D~(DK) + reaction in 
terms of the DK invariant mass is considered in Subsec¬ 
tion III ESI 


A. Elastic DK scattering amplitude 


Let us start by discussing the S -wave amplitude for 
the isospin 1 = 0 DK elastic scattering, which we denote 
T. It can be written as 


T~\s) = V-\s) - G(s) => T(s) = V(a)(l + G(s)T(s)) , 

(!) 

where G(s) is a loop function bearing the unitary or right 
hand cut, 

d 4 l 1_1_ 

7 (27t) 4 l 2 — m 2 K + ie (q — l) 2 — m 2 D + it 

and s = q 2 is the invariant mass squared of the DK 
system. This function needs to be regularized, and this 
is accomplished in this work by means of a subtraction 
constant, a(/x). In this way, the G function can be written 
as: 


m D m K 


167t 2 G(s) = a{fi) + log 

A* 

v (, s-A + v 

+ TT lo S-r-— + log 

2 s V —s + A is 


A m 2 D 

s los < 

s + A + v 


—s — A 


(3) 


A = 


L D 


~ rn?k i v — to d, vr? K ) 


where A(a;, y , z) = x 2 + y 2 + z 2 — 2 xy — 2 yz — 2 zx is the 
Kallen or triangle function. In Eq. Q, V{s) is the poten¬ 
tial, typically extracted from some effective field theory, 
although a different approach will be followed here (see 
below). 

The amplitude T(s) can also be written in terms of the 
phase shift <5(s) and/or effective range expansion param¬ 
eters, 


rru \ 

1 (s) =---:- 

PK ctgd - IPK 


8tt^ 

1 I 2 

- + T^I'OPk ^ WK 
a l 


(4) 


with 


Pk(s) 


A 1 G(s ; Mj-., Mjj) 
2 yfs 


(5) 
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the momentum of the K meson in the DK center of mass 
system. Above, a and r 0 are the scattering length and 
the effective range, respectively. 

In this channel and linked to it we find the .D*)}" (2317) 
resonance, the object of study of this paper, below the 
DK threshold, the latter being located roughly above 
2360 MeV. This means that the amplitude has a pole at 
the squared mass of this state, M 2 = so> so that, around 
the pole, 


T( S ) 



+ regular terms, 


( 6 ) 


being g the coupling of the state to the DI\ channel. 
From Eqs. (CD and ©, we see that (the following deriva¬ 
tives are meant to be calculated at s = so): 

1 dT~\s) dV~\s) dG(s) 

g 2 ds ds ds 

We have thus the following exact sum rule, 


1 


= 9 


dG(s) 

ds + ds ) ' 


( 8 ) 


In Ref. [34| it has been shown, as a generalization of the 
Weinberg compositeness condition [29| (see also Ref. 


and references therein), that the probability P of finding 
the channel under study (in this case, DI\) in the wave 
function of the bound state is given by: 


P=-g 


i dG(s) 

ds 


(9) 


In this way, the quantity ag 2 represents the probabil¬ 
ity of finding other components beyond DK in the wave 
function of (2317). The following relation can also 
be deduced from Eqs. Cl and ©: 


1 — Pdk dG(s) 
Pdk ds 


S=S 0 


(14) 


We can now link this formalism with the results of 
Ref. [HJ], where a reanalysis is done of the energy levels 
found in the lattice simulations of Ref. [27} ■ In Ref. [H| , 
the following values for the effective range parameters are 
found: 


a 0 = —1.4 ± 0.6 fm , r 0 =-0.1 ± 0.2 fm . (15) 

Also, in studying the D*o"(2317) bound state, a binding 
energy B = Md+Mk — M d *+ = 31±17 MeV is found in 
Ref. [28|. The probability Pdk is also studied, and the 
value Pdk = 0.72±0.12 is found. Hence, for our analysis, 
in which synthetic data for the reaction B® —> ( DI\) + D~ 
will be generated, we can start from the hypothesis that 
a bound state exists in the DK channel, with a mass 
M d *+ = 2317 MeV (the nominal one), and with a prob¬ 
ability Pdk = 0.75. This implies, from Eq. (ThU) . the 
value a = 2.06-10 -3 GeV~ 2 . Finally, for the subtraction 
constant in the G function, Eq. ©, we shall take, as in 
Ref. [HI, the value a(g) = —1.3 for g = 1.5 GeV. Note 
that dG(s)/ds does not depend on g or a(g). 


B. Decay amplitude and invariant DK mass 
distribution in the B° —> D~ {DK) + decay 


while the rest of the r.h.s. of Eq. © represents the prob¬ 
ability of other channels, and hence the probabilities add 
up to 1. Finally, if one has an energy independent poten¬ 
tial the second term of Eq. © vanishes, and then P = 1, 
that is, the bound state is purely given by the channel 
under consideration. In Ref. [34|, these ideas are gener¬ 
alized to the coupled channels case. 

Let us now apply these ideas to the case of DK scat¬ 
tering. From Eq. CD it can be seen that the existence of 
a pole implies 


V (s) ~ G(s 0 ) + a(s — So) + • ■ 
dV-^s) 


ds 


( 10 ) 

( 11 ) 


in the neighbourhood of the pole. If we assume that the 
energy dependence is smooth enough to allow us to retain 
only the linear term in s, we can now write the amplitude 
as 


T~ 1 (s) = G(s 0 ) — G(s) + a(s — s 0 ) , (12) 


We now discuss the amplitude for the decay B° s —>• 
D~{DK) + decay, and its relation to the DI\ elastic scat¬ 
tering amplitude studied above. The basic mechanism for 
this process is depicted in Fig. [U where, from the sb ini¬ 
tial pair constituting the B®, a cs pair and a sc pair are 
created. The first pair produces the D ~, and the DK 
state arises from the hadronization of the second pair. 
Let us consider in some more detail the hadronization 
mechanism. To construct a two meson final state, the cs 
pair has to combine with another qq pair created from 
the vacuum. We introduce the following matrix, 



/ uu ud us uc \ 

du dd ds dc 

su sd ss sc ’ 

\ cu cd cs cc j 

which fulfils: 

M 2 = (vv)(vv) = v(vv)v 

= (uu + dd + ss + cc) M . 


(16) 


and the sum rule in Eq. © becomes: 

Pdk = 1 — ag 2 ■ 


(13) 


(17) 
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The first factor in the last equality represents the qq cre¬ 
ation. This matrix M is in correspondence with the me¬ 
son matrix (f>: 


( JL 

V3 


V 


' V2 ' 

7T~ 

K~ 

D° 


v 

Ve 


v 

Vs 


7T 

77 ° 

' V2 ‘ 
K° 
D+ 


V 

V6 


K+ 

K° 


D° \ 
D~ 


VM. _2_ n- 

% 


7 

( 18 ) 

The hadronization of the cs pair proceeds then through 
the matrix element (A7 2 ) 43 , which in terms of mesons 
reads: 


(<^) 43 = K+D° + K°D+ + • • • , (19) 


where only terms containing a KD pair are retained, 
since coupled channels are not considered in this work. 
We note that this KD combination has 7 = 0, as it 
should, since it is produced from a cs, which has 7 = 0, 
and the strong interaction hadronization which conserves 
isospin (the qq with the quantum numbers of the vacuum 
has 7 = 0). 

Let t be the full amplitude for the process B ° —>■ 
D~ (DK ) + , which already takes into account the final 
state interaction of the DK pair. Also, let us denote by 
v the bare vertex for the same reaction. To relate t and 
v, that is, to take into account the final state interaction 
of the DI\ pair, as sketched in Fig. [2j we write: 


t = v + vG(s)T(s) = v(l + G(s)T(s)) . (20) 


From Eq. <m>, the previous equation can also be written 
as: 


t = v 


T(s) 
V(s) ' 


( 21 ) 


Because of the presence of the bound state below thresh¬ 
old, this process will depend strongly on s in the kinemat- 
ical window ranging from threshold to 100 MeV above it, 
so we can safely assume that t depends only on s. Hence, 
the differential width for the process under consideration 
is given by: 


dr 

dy/s 


1 

32tt 2 M|„ 


P d -Pk 


C Pd~Pk 


T(s) 

V(s) 


, ( 22 ) 


where the bare vertex v has been absorbed in C, a global 
(but otherwise not relevant) constant, and where px is 
given in Eq. ([5]) and p D - is the momentum of the Dj 
meson in the rest frame of the decaying B 3 , given by: 


A 1 / 2 (M 2 


POT = 




M 2 d -,s ) 


2 M 


bo 


(23) 


III. RESULTS 

We want to investigate the presence of the (2317) 
state in the ( DI\) + scattering amplitude. In order to 



Central Value 5 % 10 % 

10 3 a (GeV" 2 ) 

9 +0.17 +0.10 

z.uu -0.40 -1.09 

M D * n (MeV) 

9*317 

z,ot/ —24 -73 

a{p) 

on +0.15 +0.27 

i.OU -0.37 -0.49 

\9\ (GeV) 

i i n +1.0 +2.2 

ii>u -0.6 -1.1 

a 0 (fm) 

n +0.2 +0.4 

-0.2 -0.5 

ro (fm) 

n i a +0.06 +0.16 

u ' 14; -0.03 -0.04 

Pdk 

r\ yr. +0.07 +0.16 

u -'° -0.06 -0.11 


TABLE I: Fitted parameters (a, A7 n .+ and a(/d) and pre- 

s 0 

dieted quantities (|<?|, ao, ro, Pdk) for p = 1.5 GeV. The 
second column shows the central value of the fit, whereas the 
third (fourth) column presents the errors (estimated by means 
of MC simulation) when the experimental error is 5% (10%). 


explore the sensitivity of the decay B ° —»• ( DK) + D~ 
to the presence of this bound state, we generate syn¬ 
thetic data from our theory for the differential decay 
width for the process with Eqs. (1^1) and m- We gen¬ 
erate 10 synthetic points in a range of 100 MeV start¬ 
ing from threshold. To each centroid, we assign the 
value obtained with the central values explained in Sub¬ 
section III Al (10 3 a = 2.06 GeV~ 2 , a(p) = —1.3, and 
M d *+ = 2317 MeV). We shall study two different cases, 
in which each experimental point is given an error of a 
5% or a 10% of the highest value of the differential de¬ 
cay width. Taking these synthetic data as experiment- 
given data, we perform the inverse problem of analysing 
them with our theory. Obviously, the reproduction of 
the data must be perfect, but we recall that the scope 
here is to investigate the experimental accuracy that is 
actually needed to obtain reliable values for the quanti¬ 
ties fitted or predicted from our theory (M D *+, a o, ro, 
and Pdk)- The analysis of these synthetic data goes as 
follows. We generate around 2-10 3 sets of random exper¬ 
imental points, in which each centroid is varied around 
its theoretical value according to a Gaussian distribution 
with the error given to each point. For each of these sets 
of random points, the parameters are fitted to the data. 
After the whole run, a central range, containing a 68% 
of the values of the considered quantities (the differential 
decay width, the fitted parameters, and the predicted 
values) is retained. It is worth stressing here that, since 
the centroid of the experimental point in each set of ran¬ 
dom experimental points is varied, a good reproduction 
of the random synthetic data is quite natural but not 
completely trivial. 

The generated synthetic data are shown in Fig. [3| As 
explained, they have two different error bars, the smaller 
one corresponding to a 5% experimental error and the 
larger one to a 10%. As commented above, they exactly 
match the central curve (dash-dotted line) produced with 
the central parameters of the theory. A solely phase space 
distribution (i.e., a differential decay width proportional 
to p d ~pk, but with no other kinematical dependence 
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FIG. 2: Diagrammatical interpretation of Eq. @DJ, in which DK final state interaction is taken into account for the decay 
B° —► D^{DK) + . The dark square represents the amplitude t for the process, in which the final state interaction is already 
taken into account. The light square represents the bare vertex for the process, denoted by v. Finally, the circle represents the 
hadronic amplitude for the elastic DK scattering. 



FIG. 3: Differential decay width for the reaction B® —> 
DJ (DK) + . The synthetic data (generated as explained in the 
text) are shown with black points. The smaller (larger) error 
bars correspond to a 5% (10%) experimental error. The dash- 
dotted line represents the theoretical prediction obtained with 
the central values of the fit. The light (dark) bands corre¬ 
spond to the estimation of the error (by means of a MC sim¬ 
ulation) when fitting the data with 5% (10%) experimental 
error. The dashed line corresponds to a phase space distribu¬ 
tion normalized to the same area in the range examined. 


of dynamical origin) is also shown in the figure (dashed 
line). The first important information to be extracted 
from the figure is that the data are clearly incompatible 
with this phase space distribution. This points to the 
presence of a resonant or bound state or, at least, to 
some strong final state interactions. Two error bands 
are shown in the same figure, the lighter and smaller 
(darker and larger) one corresponding to a 5% (10%) 
experimental error. The fitted parameters (a(/z), M D *+, 

and a) are shown in Tabled together with their errors. 2 
Note that, with a 5% experimental error, we get M D ,+ = 


2 To avoid unphysical values of the fitted parameters a(fi) and 
a, which could numerically reproduce each set of the randomly 
generated experimental points, they are restricted to vary within 
a sensible range, but making sure that this range is larger than 
the error obtained for these two parameters and shown in Table[I] 


2317124 MeV, and if the error is increased to 10%, the 
value is M D *+ = 2317 I 73 MeV. Here we are mainly 
concerned with the upper error, in the sense that it is the 
one that defines if the bound state is clearly below the 
DK threshold (which is slightly above 2360 MeV) or not. 
Considering this error, we see that the mass obtained is 
well below the threshold, at the level of 2cr (3<r) for the 
case of a 5% (10%) experimental error. This is a good 
information: experimental data with a 10% error, which 
is clearly feasible with nowadays experimental facilities, 
can clearly determine the presence of this below threshold 
state D*q“(2317). 

We can also determine Pdk , the probability of finding 
the DK channel in the -0*^(2317) wave function. It is 
shown in the last row of Table |U As stated, the central 
value Pdk = 0.75 is the same as the initial one, but we 
are here interested in the errors, which are small enough 
even in the case of a 10% experimental error. This means 
that with the analysis of such an experiment one could 
address with enough accuracy the question of the molec¬ 
ular nature of the state (D*o"(2317), in this case). 

Finally, it is also possible to determine other parame¬ 
ters related with DK scattering, such as the scattering 
length (ao) and the effective range (ro). They are also 
shown in Table U They are compatible with the lattice 
QCD studies presented in Refs. j27j, [28J. Namely, the 
results from Ref. [28| are shown in Eqs. m, and their 
mutual compatibility is clear. 


IV. CONCLUSIONS 

In the present work we have selected a reaction which 
is both Cabbibo and color favored, the B° s —> Dj(DK) + 
weak decay, and have looked at the DK invariant mass 
distribution from where we expect to obtain relevant 
information on the nature of the U*q“( 2317) resonance 
when actual data are available. For this purpose we have 
performed a simulation of the experiment taking informa¬ 
tion from experiment about the mass of the Z?*q“( 2317) 
resonance and from a recent QCD lattice analysis on 
an analytical representation of the KD scattering arnpli- 
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tude. This information has served us to make predictions 
on the shape of the KD invariant mass distribution close 
to the KD threshold. After that we have taken these 
results and we have assumed they are actual ’’experimen¬ 
tal data”, associating to them an ’’experimental error” of 
5% or 10%. Then we have made a fit to these ’’synthetic 
data” in order to extract from there the KD scattering 
amplitude, above and below threshold. We prove that 
with both errors, typical of present experimental data of 
spectra in B decays, one can obtain the KD scattering 
amplitude with enough precision to predict that there 
is a KD bound state. We also predict the scattering 
lenght and effective range of the KD interaction and, 
very important, we show that we can predict, with rel¬ 
atively small error, the probability of the mesonic KD 
component in the wave function of the J5 *q"( 2317) reso¬ 
nance. From the QCD lattice results one induces about 
70% probability and we show that this number can be 
obtained from the analysis of the B decay spectra with 
sufficient precission to make the number significative of 
the main nature of the -0*^(2317) resonance as a basi¬ 
cally I\D molecular state with a smaller mixture of other 
components. 

The study done here should stimulate the implementa¬ 


tion of the experiment, for which we have made estimates 
of a relatively large branching fraction. 
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